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Abstract 

We analyze different examples to show that the so-called generalized Hartman effect is an erroneous presumption. 
The results obtained for electron tunneling and transmission of electromagnetic waves through superlattices and 
Bragg gratings show clearly the resonant character of the phase time behavior so that a generalized Hartman effect is 
not expected to occur. A reinterpretation of the experimental results in double Bragg gratings is proposed. 
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Background 

The Hartman [1] effect is known as the independence of 
the tunneling time on the barrier width as this parame- 
ter gets large. It has been shown that the experimental 
evidences of this effect on the transmission times of pho- 
tons and electromagnetic pulses [2-6] are compatible with 
phase time calculations [7]. The Hartman effect has been 
investigated in various ways by extending the system not 
only for a single barrier but also for double [8,9] and 
multiple barrier [10,11] structures. Olkhovsky, Recami, 
and Salesi came out with the idea that for non-resonant 
tunneling through two potential barriers, the tunneling 
time (which is a phase time) is independent not only of 
the barrier width but also of the barrier separation [8]. 
The approximations introduced in this reference to obtain 
the unknown coefficients, led these authors to unphysical 
results like the generalized Hartman effect. This has been 
called the generalized Hartman effect (GHE). The two- 
barrier problem can be solved without approximations, 
see for example, in the work of Pereyra [12]. An experi- 
ment to check this effect was performed by Longhi et al. 
[10] where optical pulses of 1, 550 nm wavelength were 
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transmitted through a double-barrier system of Bragg 
gratings. In this reference, non-conclusive and inappro- 
priately presented results for five different values of the 
gratings separation were reported. Most of the theo- 
retical conclusions were based on questionable formu- 
las and unnecessarily involved calculations. For example, 
Equation 2 (used in Equations 3 and 4) of [8] is not the 
actual transmission coefficient through a double Bragg 
grating. A criticism on the mathematical rigor on GHE 
is also given by S. Kudaka and S. Matsumoto [13,14]. It 
is easy to check from a straightforward calculation, or 
from the precise and general formulas published in [7] 
as quoted below, that the phase time for a double barrier 
(DB) with separation L has the structure 



21 dk 

x = T 2 — + T 2 Ai(FsmkL + Gcos/cL) 

T dco 

+ T 2 A r (F cos kL + Gsin/cl) 



(1) 



Spri 



with T 2 and T the double- and single-barrier transmis- 
sion coefficients, respectively, k the wave number, go the 
frequency and Au A r , F, and G simple functions of the 
potential parameters (P. Pereyra and H. P. Simanjuntak, 
unpublished work). Despite this clear dependence on L, 
involved and contradictory arguments lead to establish 
that r becomes independent of L [8,10,11]. In the follow- 
ing we will consistently use a for the separation between 
barriers. 
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For the inference of a generalized Hartman effect to be 
meaningful for multi-barriers, double superlattices (SLs) 
or double Bragg gratings (BG), one would of course need 
to keep the physical parameters [like the energy (wave- 
length) of the particle (wave)] fixed as the length between 
barriers is increased. The tunneling and transmission 
times behavior should be taken with care when one tries 
to find a Hartman effect due to barrier separation in multi- 
barrier systems [8,11] since, in general, the density of res- 
onance energies grows rapidly as the separation increases. 
It is well known that the non-resonant gaps in the band 
structure of a SL or a BG become resonating when these 
systems are divided and separated; and the separation is 
increasingly varied. This was already recognized in [15] 
(for double SL) and in [10] (for double BG). On the other 
hand, it is well known that the tunneling time follows the 
resonant band structure [7,16]. Thus, it is not possible 
to keep increasing the separation between barriers and 
superlattices without crossing resonances. For this rea- 
son, visualized here with specific examples for electrons 
and electromagnetic waves, the existence of a generalized 
Hartman effect is a rather questionable issue. For these 
examples we perform first principle calculations using the 
actual transmission coefficient of the system (such as that 
of double BG in the experiment in [10]) so that we can 
justify completely that the so-called generalized Hartman 
effect is erroneous. 

To study the Hartman effect and to criticize the pre- 
sumption of a generalized Hartman effect in superlattices, 
Bragg gratings, and multi-barrier systems, we will use the 
theory of finite periodic system that allows straightfor- 
ward calculation of the phase time. For electron tunneling, 
we shall assume periodic and sectionally constant poten- 
tials with cells of length i c = a + b and a barrier of 
width b and strength V 0 in the middle. For electromag- 
netic waves, each cell consisting of dielectrics 1 and 2 will 
contain a dielectric 2 of length b in the middle. In this case 
€i, Hi, and fii (with i = 1, 2) are the corresponding permit- 
tivities, refractive indices, and permeabilities; the regions 
outside the SL are assumed to be air. For Bragg gratings, 
the refractive indices are periodic. 

Methods 

If we have a Gaussian wave packet (of electrons or elec- 
tromagnetic waves) through a SL of length ni c — a, the 
centroid phase time (which is taken here as the tunneling 
or transmission time) is given by [7,17,18] 



dk 

- aH dE + 



1 daj 

-U 2n -,(a s) — 



ai r a R "I dan I 



Here a = + iotj is the (1,1) element of the single-cell 
transfer matrix M; U n (aR) are the Chebyshev polynomi- 
als of the second kind evaluated at cxr; and a n is the (1,1) 
element of the w-cell transfer matrix M n . This is given by 
[16] 



a n = U n (a R ) -a*U n -i(a R ). 
At resonance, where U n -\ = U<m-\ = 0, we have [16] 

/ V7T > 



a R = COS 



V7T \ 

n / 



v = 0,±1,... 



(3) 



(4) 



The expression for the tunneling or transmission time 
simplifies as 



nrwtj daR 



(5) 



The tunneling time in Equation 2 is exact and general 
and valid for arbitrary number of cells, barrier width, 
and barrier separation. Thus, one can check the exis- 
tence or not of a (generalized) Hartman effect at will. 
For concrete examples, we consider superlattices like 
(GaAs/Alo.3Gao.7As) w /GaAs, with electron effective mass 
m\ = 0.067 m in GaAs layers, = 0.1 m in Alo.3Gao.7As 
layers (m is the bare electron mass) and V Q = 0.23 eV, and 
Bragg gratings with periodic refractive index. 

Results and discussion 

Electron tunneling 

If we consider electrons through superlattices with unit 
cell length i c = a + b, we will have 

a = M\\ = exp[ika] (cosh qb — i(q 2 — k 2 ) sinhqb/ 2qk) 

(6) 



with k = ^2m A E/h and q = ^2m B (V 0 -E)/h. When 
mA, rriB and V Q are taken as fixed parameters, we choose 
a = 100 A and b = 30 A. 

For a single barrier, n = 1, the tunneling time ri plot- 
ted in Figure 1 as a function of the reduced barrier width 
b/X shows the well-known Hartman effect. The energy 
E is kept fixed and X = 27th/ '^/2wiaE is the de Broglie 
wavelength. 

With aR and aj growing exponentially with the barrier 
width b, one can easily show from Equation 2 that for large 
b, the non-resonant tunneling time approaches that for a 
single barrier, i.e., x n (E) & Ti(E) as 



r n (E) 



Vo^rn A /m B 



^(V 0 -E)E V 0 - E(l - m A /m B ) 



(7) 



This is the well-known Hartman effect. Since this quan- 
tity becomes also independent of the barrier separation 
[8,11] a, it has been taken as the analytical evidence of a 
generalized Hartman effect. However, such an approxima- 
tion that leads to the independence on a and n is obtained 
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Figure 1 The single-barrier transmission coefficient and the 
tunneling time. The single-barrier transmission coefficient 1 /\a\ 2 
(gray lines) and the tunneling time x\ (dark lines) as functions of the 
reduced barrier width b/k, when the electron energies are 
E = 0.122516 eV,^ = 0.15 eVand E = 0.2 eV. In the tunneling time 
curves, the Hartman effect is evident. 



by taking the limit of large b first that is strictly speak- 
ing infinite, which makes the first barrier the only one 
that matters for the incoming wave to penetrate while the 
rest of the SL is immaterial This was also pointed out by 
Winful [9]. However, Winful [9] used an approximation: 
The transmission of the double square barrier potential to 
model the transmission through the double BG. Here, we 
present calculations using the actual transmission coef- 
ficient through the double BG. As mentioned before, for 
the generalized Hartman effect to be meaningful, it should 
not matter whatever limit we take first whether on a, b, or 
n. It turns out that a non-resonant energy region becomes 
resonant as the separation a increases (see the discussion 
on the double Bragg gratings in section 'Hartman effect in 
two Bragg gratings systems'). 

The situation is completely different for resonant tun- 
neling through a SL with large but finite barrier width b 
where Equation 5 shows that the tunneling time becomes 
x n (E) oc be 2qb (since oir and aj behave as e qb for large b). 
Thus, relatively small barrier width would be needed to 
study the effect of the barrier separation and the number 
of barriers on the tunneling time. The tunneling time for 
a relatively small barrier width is shown in Figure 2 for an 
electron (with energy E = 0.15 eV) through SLs which 
number of cells are n = 3, 4, and 6. 

Looking at oir and ai, that are oscillating functions in a, 
it is clear that it is not possible to have the tunneling time 
to be independent of the barrier separation a, by keeping 
the barrier width and number of cells fixed. Therefore, the 
so-called generalized Hartman effect is at least dubious. 
The tunneling time behavior that will be found below for 
the double BG is easy to understand here. Starting with 
a certain barrier separation a, a non-resonant phase time 
becomes a resonant one as a is increased, while the other 
parameters are kept fixed. This is shown in Figure 3 where 
the tunneling time is plotted as a function of the reduced 
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Figure 2 The tunneling time r n as a function of the reduced 


barrier width. The tunneling time r n as a function of the reduced 


barrier width b/k for electrons (with energy E = 0.1 5 eV) through 


superlattices with n = 3,4, and 6. 





barrier separation, a/k, for fixed b, n, and electron energy 
E. This result shows that in this kind of systems, the 
presumption of a generalized Hartman effect is incorrect. 

The Hartman effect as a consequence of varying the 
number of cells was already discussed in [7] . In Figure 4 we 
show three qualitatively different examples on the behav- 
ior of the tunneling time as a function of n. In Figure 4a 
for energies in the gap (E = 0.15 eV and E = 0.2 eV), the 
saturation of the tunneling time exhibits the well-known 
Hartman effect. In Figure 4b, the energy lies at the edge of 
a resonant region. The phase time x n resonates for multi- 
ples of n = 21. This behavior is clearly understood if we 
consider Equations 4 and 5. Equation 4 implies that the 
same resonance energy E^ is found for different number of 
cells as long as the ratio v/n is constant. This means that 

E\ = E^ n = E\ n = From Equation 5, it is also evident 

the linear dependence of x n on n. 

The Hartman effect and the electromagnetic waves 

Electromagnetic waves have been used for discussions on 
the Hartman effect [9]. For a superlattice L(H/L) n made 
of alternating layers with refractive indices ni and nn> 
the phase time (PT) for each frequency component of a 
Gaussian wave packet through a SL of length ni c — a is 
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Figure 3 The tunneling time t& as a function of reduced barrier 
separation and fixed barrier width. The tunneling time r 6 as a 
function of reduced barrier separation a/k for fixed barrier width b, 
number of cells n = 6 and electron energy E = 0.15 eV with the 
corresponding de Broglie wavelength k. 
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Figure 4 The tunneling time x n as the number of cells n in a SL is 
varied, (a) Saturation of r n for electron energies E = 0.1 5 eV and 

E = 0.2 eV in the gap. (b) The energy is close to a resonant band-edge. 
In this case, more resonances appear as n is increased with the energy 
fixed. No Hartman effect can be inferred from this figure. 



also obtained from Equation 2 with I<l,h = ^l,hI c an d 
with [7] 



an = cos kia cos kub 



ni + n 



2n H n L 



— sin kia sin kub, (8) 



Oil 



yp- _|_ y] 2 

cos kid sin kub + — sin kia cos kub. (9) 



2n H n L 



To see the effect of varying the size of the SL on the PT, 
one has to be sure that such variation will still keep the 
wavelength inside a photonic band gap. It was shown that 
by increasing the number of cells, for fixed thicknesses of 
layers and wavelength in a gap, the PT exhibits [7] the 
observed Hartman effect [2,3]. However, this condition 
will not be possible by varying arbitrarily the thicknesses 
of the layers. The reason is that there is only a small range 
of thicknesses that one can use to keep the chosen wave- 
length to lie in a gap before going out of it and may even 
reach resonances, as shown in Figure 5 where the PT oscil- 
lates (with a band structure) and grows as a function of the 
reduced thicknesses a/k and b/k. This is analogous to the 
electron tunneling time shown in Figure 3. 

Hartman effect in two Bragg gratings systems 

We now consider the system that was taken in [10] as 
thought to support the idea of a generalized Hartman 
effect: the double Bragg gratings (DBG). Independent of 
the approximate method used in that paper, we find that 
assuming sin(/<B<2) = 0 (the only way to obtain the 
reduced expressions of Table 1 in [10]) and still consid- 
ering a as a variable are incongruous. Moreover, the idea 
that the PT becomes independent of a is incompatible 
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Figure 5 The phase times r n as functions of the reduced 
thicknesses. The phase times r n as functions of the reduced 
thicknesses o/X (a) and b/X (b) for a superlattice (L/H/L) n made of 
fused silica (L) (ni = 1 .41 ) and titanium oxide (H) (n^ = 2.22) for 
X = 702 nm.The oscillating behavior of r n with clear band structure is 
far away of any Hartman effect. 



with the Equation (4b) in their work, where a linear depen- 
dence on a is reported. In the DBG, the gratings of length 
L 0 and refractive index n(z) = no + n\ cos(2/cbz) are sep- 
arated by a distance a. The values of a considered in the 
experiment are indicated by arrows in Figure 6. The BG 
wave equation 



dz 2 



+ k 2 n 2 (z)Mz)=0> 



(10) 



when ignoring the (ni/no) 2 term for ni/no 1 
(as in [10]), becomes the Mathieu equation, in which 
solutions tyi{z) = Se(u,v;I<BZ + tt/2) and 1^2 (z) = 
So(u, v;I<bZ + 71 12) are Mathieu functions [19] with u = 
(1 + 2^i /no) Wq/c 2 //eg and v = 2 /s fnon\k/k^. The real and 
imaginary parts of the (1,1) element of the transfer matrix 
are 



a R 
a/ 



cos kriQa sin kn^a 



W 2 
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sin kriQa 
W 2 



(11) 
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Figure 6 The phase time as a function of the Bragg gratings separation, (a) The phase time as a function of the separation a between two 
Bragg gratings, for incident A. = 1 , 542 nm, k% = 6.1 074//im, no = 1 .452, ni /no = 1.8x1 0~ 4 , and L 0 = 8.5 mm. (b, c) The PT is plotted as a 
function of ca, for a = 42 mm. The phase time in (b) is the same as that in (c) but plotted from 0 to 10 ns to compare with Figure 2 in [10]. Arrows 
indicate the as in [10]. 



with W the Wronskian and Here 0\ = 9(L o ,0), 9 2 = 0(2L o + a,L 0 + a) analogously 

for Xi,2. Mi,2. vi j2 , with 

A = e 2 0\ + mvi + V2M1 + X2X1 0(z 2 ,zi) = fifa) — -j- V2KZ2) — ^ — 

B = lk ( ° 2Vl + V2Xl) ~ knMx + X2Ml) xta,*i) = - ^ dJ ^ 

C = knoifofii + M2X1) - y~ (V2O1 + X2V1) Mfe^i) = ^1(21)^2(32) - ^2(^1)^1(^2) 

22 1 v \Z2> z l) = ♦ 

= 6> 2 Xi + X2#i - k n 0 /x 2 /xi - — — ; r V2V1. (12) dz dz 

* 2 *o (13) 
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Using parameters of Longhi et al. [10] for no, n±, I<b, and 
L 0 , the non-resonant gap becomes resonant as the gratings 
separation increases. Though details are beyond the pur- 
pose of this paper, we plot in Figure 6 the PT as a function 
of the separation a for incident-field wavelength k = 1542 
nm, and as a function of the frequency co, for a = 42 mm. 
Recall that in [10], k — 1, 550 nm was considered. While 
the PT appears completely in graph (c), in (b) it is plotted 
in a different range to compare with the experiment. The 
resonant behavior of the PT with a and the absence of any 
generalized Hartman effect are evident. Similar results are 
obtained when A. = 27r/£g. 

Conclusion 

We have shown that the presumption of generalized Hart- 
man effect for tunneling of particles and transmission of 
electromagnetic waves is not correct. 
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